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1.  INTRODUCTION 


Compton  electrons  generated  in  the  air  by 
gamma  rays  produced  by  nuclear  explosions  are 
known  to  cause  transient  electromagnetic  fields 
called  the  electromagnetic  pulse  (EMP)  or  radio¬ 
flash.  The  average  forward  range  of  the  Compton 
electrons  in  the  air  is  a  principal  factor  governing 
the  strength  of  the  EMP  caused  by  bursts  in  the 
lower  atmosphere  within  about  10  km  of  the  earth’s 
surface.  For  bursts  at  high  altitudes,  a  different 
treatment  of  the  trajectory  of  the  Compton  elec¬ 
trons  is  necessary. 

For  low- altitude  bursts,  previous  estimates1  of 
the  average  forward  range  of  Compton  electrons 
produced  by  collimated  monochromatic  gamma 
rays  were  based  on  the  use  of  an  energy-indepen- 
dent  electron  transmission  factor  for  the  Compton 
electrons.  The  transmission  factor  represents  the 
cumulative  effect  of  small-  and  large-angle  scatter¬ 
ing  of  the  electrons,  causing  the  electrons  to 
straggle.  Use  of  an  energy- independent  transmis¬ 
sion  factor  neglects  both  the  true  average  Compton 
recoil  energy  and  the  recoil  energy  dependence 
upon  the  scattering  angle. 

In  this  investigation,  the  average  forward 
range  is  computed  using  an  energy-dependent 
transmission  factor  in  connection  with  the  Comp¬ 
ton  angular  scattering  distribution.  This  is  done  for 
a  range  of  gamma  ray  energies  from  0.3  to  10  MeV. 
Simple  numerical  fits  to  the  results  are  included. 


2.  METHOD  OF  CALCULATION 

The  well-known  Klein-Nishina  differential 
scattering  cross  section  for  Compton  scattering  is 

A(ae)  =  2*  .in  8  d B  1+a(1l_C08^  ] * 

(l+COB1  $\  (  ,  0^(1  —  CO*  8)1  \ 

2  /  V  d+co#1  0)(l+a<l  —cos  0))/ 

per  electron, 

John  S.  Mmtik,  7V  Compton  Current  Lot  Alamo*  Scientific  Laboratory, 
Defense  Atomic  Support  Agency,  1882-1  (5  November  186$). 


where 

6  =  gamma  scattering  angle, 
or  =  gamma  energy  in  units  of  0.51 1 
MeV,  and 

r0  =  classical  electron  radius. 

The  electron  recoil  angle  0  is  given  by 
cot  0  —  (1+a)  tan  (6/2)  . 

The  scattered  gamma  energy  E  is 
Ey  =  Eyll+a(l  —  cos  0))  , 

where  Ey  =  0.511a  is  the  incident  gamma  en¬ 
ergy  (in  MeV).  The  electron  recoil  energy  E  is 

E  =  Ey— Ey  . 

The  range  of  the  Compton  electron  is  commonly 
given  in  two  alternative  forms:  (1)  the  practical  or 
extrapolated  range  or  (2)  the  maximum  or  CS DA 
range  (i.e.,  continuous  slowing  down  approxima¬ 
tion).  The  former  is  obtained  from  experiments 
and  the  latter  from  slowing-down  theory.  The 
extrapolated  range  is  commonly  estimated  to  be 
about  0.8  of  the  maximum  range.  The  difference 
is  attributed  to  multiple  scattering  of  the  elec¬ 
trons.  Results  were  obtained  in  this  investigation 
for  an  example  of  each  form  for  electrons  in  air 
of  density  p : 

( 1 )  pR**  =  0.368 EA,  A  =  1 .26-0.064  InE  (cm) , 

(2)  pR*  *  0.6 1  E7(E+0.26)  (cm) , 

where  the  first  formula  is  a  frequently  used 
approximation  and  where  the  latter  formula  is  due 
to  Longmire  as  cited  by  Malik.1  As  suggested. 

Re*  *  0.8  R^  . 

The  electron  transmission  factor  (also  called 
“transmitted  number  fraction”)  has  been  mea¬ 
sured  by  Ebert  et  al,2  who  suggested  the  following 
expression  for  the  transmission  factor  T(x) : 

'P  J.  Ebert  A.  F.  Lemon,  ond  E  V  Lent  Lawrence  Rodiotton  Loborotory , 
UCRL-71462  ( December  1968 %  o*  cited  by  J.  A.  Lonergan  ond  D.  C  Shreve. 
Parametric  Fit  to  Electron  Trmntport  Prapertic*.  Proeeeduqt  of  the  Notional 
Symposium  on  Natural  mod  Monmode  Radiation  in  Spore,  E  A.  Worman,  rd. 
National  Aeronautics  ond  Space  Administration,  TM  X-2440  (January  1972). 
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TABLE  1.  AVERAGE  ELECTRON  RANGE 
R  AS  A  FRACTION  OF  EXTRAPOLATED 
RANGE  R„,  USING  EBERT'S  FORMULA 


T(x)  =  exp  ( —  A(x/Rejr)B]  , 
where 


A  =  (1  —  1/B)1  -B  , 

B  =  f  387  E  1  1/4  .  and 

l  Z(1+7.5X10~5  ZE:)  J 

Electron  energy  E 
(MeV) 

Ratio  R/R*.* 

0.03 

0.7488 

0.1 

0.6243 

Z  =  the  atomic  number  (7.2  for  air). 

0.5 

0.6404 

1.0 

0.6650 

The  following  integrals  define  Compton  total 

1.5 

0.6812 

cross  section  ( ae ),  average  Compton  recoil  energy 

2.0 

0.6932 

(E),  average  forward  direction  cosine  of  the  Comp¬ 

2.5 

0.7027 

ton  electron  (cos  0),  and  Compton  electron  aver¬ 

3.0 

0.7104 

age  forward  range  (R): 

3.5 

0.7170 

5.0 

0.7320 

ro=7T 

7.0 

0.7458 

*  ~)«-o 

10.0 

0.7597 

E 


COS  0 


R 


i 

a »  3  0=0 

if'- 

ae  3  0=0 

L  {e=n 

O'  3  0=o 


E  d(ae)  , 


cos  <f>  d(aP)  , 


ft(E)  COS  <t>6(ae)  . 


The  average  electron  range  R  is 


ft(E)  —  R** 


dTW 

dx 


dx  , 


independent  values  of  2/3  and  0.63  have  been 
used  previously  as  cited  by  Malik.1 


It  will  be  convenient  to  express  the  Compton 
electron  average  forward  range  R  for  a  gamma 
energy  Ey  as 

R  =  C(Ey)  Re*  ( E )  cos0  * 

where  the  extrapolated  range  R*,*  is  determined 
from  the  average  Compton  electron  energy  E ,  the 
Compton  electron  average  forward  direction  cosine 
cos  0  is  determined  from  the  gamma  energy,  and 
the  factor  C  is  determined  from  the  gamma  energy. 
Once  R,  E,  and  cos  0  are  known,  f  is  readily 
found. 


which,  using  the  Ebert  form  for  T(x) ,  is  easily 
found  to  be 

ft(E)  -  ReTaar(2— a)  , 
where 

a  -  1  —  1/B  . 

The  ratio  ft/Re*  is  tabulated  for  a  few  values 
of  electron  energy  E  in  table  1.  Energy* 


3.  RESULTS  OF  CALCULATION 


The  integrals  for  ae  ,  E ,  cos  0 ,  and  R  were 
integrated  numerically  for  a  range  of  gamma  ray 
energies  from  0.01  to  10.0  MeV.  This  was  done 
both  for  the  extrapolated  range  R**  and  the  esti- 

John  5.  Molik.  The  Compton  Current.  Los  Ahmos  Scientific  Lmbormtor\ 
Drfmsr  Atomic  Support  A/trncy.  1882-  fS  \ov*mbrr  /965>. 
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mated  extrapolated  range  0.81^,,  since  expressions 
for  and  0.8  yield  somewhat  different  values. 
Results  are  shown  in  table  2  for  R**  and  in  table  3 
for  0.8  Rq.  In  addition  to  E,  cos0,  and  R,  theresults 
include  the  factor  £,  and  also  the  value  for  R  that  is 
obtained  if  an  energy-independent  ratio 
ft/R**  =  2/3  is  used. 

Comparison  of  R  with  R(R/Re*  =  2/3) 
shows  the  effect  of  neglecting  the  energy  depen¬ 
dence  of  the  transmission  factor.  The  difference  is 
negligible  for  a  gamma  energy  of  1.5  MeV,  and 
increases  the  more  the  energy  level  differs  from  1 .5 
MeV.  For  example,  the  difference  is  7  percent  for 
gamma  energies  of  0.3  MeV,  5  percent  for  0.7 


MeV,  5  percent  for  3^5  MeV,  and  10  percent  for  7 
MeV.  The  value  for  R(R/Rejr  =  2/3)  is  too  low  at 
higher  gamma  energies  and  at  very  low  gamma 
energies.  This  is  understandable  since  ( 1 )  at  high 
gamma  energies  most  recoil  electrons  are  scattered 
nearly  forward  with  energies  so  high  that 
ft/R«  =  2/3  is  an  underestimate,  and  (2)  at 
very  low  gamma  energies  for  which  Compton 
recoil  energies  approach  0.02  MeV,  the  Ebert 
formula  for  the  transmission  factor  begins  to  fail 
(predicting  an  increase  in  R/R**  to  a  value  of  0.749 
at  E  =  0.03  MeV,  for  example).  The  failure  in 
the  Ebert  formula  is  due  to  the  quantity  1/B 
approaching  and  exceeding  unity.  For  this  reason, 
the  computation  cannot  reasonably  be  applied  to 
gamma  energies  less  than  about  0.3  MeV. 


TABLE  2.  COMPUTED  RESULTS  USING  ELECTRON  EXTRAPOLATED  RANGE" 


Ey 

(MeV) 

E 

(MeV) 

pR** 

(g/cm2) 

COS  0 

pR 

(g/cm2) 

pRfor  ft/R,,  =  2/3 
(g/cm2) 

0.01 

2.181 

0.188(— 3) 

0.666(— 7) 

0.656 

0.953(— 7) 

0.635(-7) 

0.1 

1.461 

0.0138 

0.516(— 3) 

0.654 

0.492(— 3) 

0.383(— 3) 

0.15 

1.137 

0.0272 

0.171(— 2) 

0.658 

0.128(-2) 

0.122(— 2) 

0.2 

1.025 

0.0433 

0.374(— 2) 

0.663 

0.255(-2) 

0.261(-2) 

0.3 

0.937 

0.0809 

0.0103 

0.677 

0.655(— 2) 

0.699(— 2) 

0.4 

0.897 

0.124 

0.0200 

0.692 

0.0124 

0.0133 

0.5 

0.873 

0.171 

0.0325 

0.707 

0.0200 

0.0214 

0.6 

0.857 

0.221 

0.0473 

0.721 

0.0292 

0.0309 

0.7 

0.844 

0.273 

0.0643 

0.733 

0.0398 

0.0418 

0.8 

0.834 

0.327 

0.0831 

0.745 

0.0516 

0.0539 

1.0 

0.820 

0.440 

0.125 

0.766 

0.0787 

0.0811 

1.5 

0.799 

0.742 

0.251 

0.807 

0.162 

0.162 

2.0 

0.787 

1.062 

0.397 

0.835 

0.261 

0.257 

2.5 

0.781 

1.393 

0.555 

0.856 

0.371 

0.360 

3.0 

0.777 

1.732 

0.721 

0.872 

0.489 

0.468 

3.5 

0.775 

2.077 

0.893 

0.885 

0.612 

0.580 

5.0 

0.771 

3.140 

1.431 

0.911 

1.005 

0.931 

7.0 

0.771 

4.597 

2.167 

0.931 

1.556 

1.412 

10.0 

0.773 

6.836 

3.273 

0.948 

2.399 

2.133 

aSymbols  are  defined  as  follows, 

Ey  *=  gamma  energy 

*  correction  factor  obtained  for  pR  =*  Ce*pRe*(E)  cos  0 
E  =  average  recoil  electron  energy 
pR^  *  extrapolated  range  for  electron  energy  E,  air  density  p 
cos  0  —  Compton  electron  forward  direction  cosine 

pR  —  Compton  electron  average  forward  range  for  transmission  factor  T(E),  air  density  p 
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It  is  emphasized  that  these  results  for  the 
average  forward  range  of  Compton  electrons  are 
valid  only  for  incident  collimated  monochromatic 
gamma  rays  undergoing  a  first  scatter.  Subsequent 
scatters  would  occur  at  lower  gamma  energies  and 
in  different  directions  than  for  the  first  scatter. 
However,  accurate  Compton  electron  current  scores 
could  be  obtained  from  multiple-scatter  gamma 
transport  calculations  if  the  Compton  electron 
average  forward  range  presented  here  were  used  to 
score  each  Compton  scatter. 

The  expression 

R  =  C(Ey)R«*(E)co3<A 


uses  the  fit  to  the  extrapolated  range  R^fE )  in  air 
for  the  average  Compton  recoil  energy  E ;  in  this 
expression,  the  factor  £(Ey),  a  function  of  the 
incident  gamma  energy  Ey  (in  MeV),  is  approx¬ 
imated  to  within  ±0.003  on  the  interval  0.3  MeV  < 
Ey  <  10  MeV  by 

CeX  -  1.22/(l  +  1.2065x-0.626x2)  , 

where 

x  =  Ey/(Ey+0.7)  . 

When  the  extrapolated  range  is  estimated  as 
0.8Ro(Ey)  where  is  Longmire’s  fit  to  the  CSDA 
range,  tne  factor  £  is  approximated  to  within 
±0.003  on  the  same  interval  by 


TABLE  3.  COMPUTED  RESULTS  USING  ELECTRON  CSDA  RANGE" 


Ey 

(MeV) 

Co 

E 

(MeV) 

0.8pRo 

(g/cm2) 

COS  0 

pR 

(g/cm2) 

pRfor  R/R**  =  2/3 
(g/cm2) 

0.01 

1.877 

0.188(— 3) 

0.660(— 7) 

0.656 

0.813(-7) 

0.542(-7) 

0.1 

1.569 

0.0138 

0.339(— 3) 

0.654 

0.348(— 3} 

0.272(— 3) 

0.15 

1.234 

0.0272 

0.126(— 2) 

0.658 

0.102(— 2) 

0.984(— 3) 

0.2 

1.112 

0.0433 

0.301(— 2) 

0.663 

0.222(— 2) 

0.229(-2) 

0.3 

1.000 

0.0809 

0.937(— 2) 

0.677 

0.635(— 2) 

0.678(-2) 

0.4 

0.940 

0.124 

0.0195 

0.692 

0.0127 

0.0136 

0.5 

0.900 

0.171 

0.330 

0.707 

0.0210 

0.0224 

0.6 

0.872 

0.221 

0.494 

0.721 

0.0310 

0.0328 

0.7 

0.851 

0.273 

0.0681 

0.733 

0.0425 

0.0447 

0.8 

0.835 

0.327 

0.0889 

0.745 

0.0553 

0.0577 

1.0 

0.813 

0.440 

0.135 

0.766 

0.0841 

0.0866 

1.5 

0.786 

0.742 

0.268 

0.807 

0.170 

0.171 

2.0 

0.775 

1.062 

0.416 

0.835 

0.269 

0.265 

2.5 

0.770 

1.393 

0.573 

0.856 

0.378 

0.366 

3.0 

0.768 

1.732 

0.735 

0.872 

0.492 

0.471 

3.5 

0.767 

2.077 

0.901 

0.885 

0.612 

0.580 

5.0 

0.769 

3.140 

1.415 

0.911 

0.991 

0.918 

7.0 

0.773 

4.597 

2.123 

0.931 

1.528 

1.387 

10.0 

0.779 

6.836 

3.214 

0.948 

2.375 

2.111 

aSymbols  are  defined  as  follows . 


Ey  =  gamma  energy 

=  correction  factor  obtained  for  0.8pR  =  ^o0.8pRo(E)  cos 
E  =  average  recoil  electron  energy 

0.8pRo=  estimated  extrapolated  range  from  CSDA  range  Rq  for  electron  energy  E  9  air  density  p 
cos  0  —  Compton  electron  average  forward  direction  cosine 

pR  =  Compton  electron  average  forward  range  for  transmission  factor  T(E),  air  density  p 
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Co  *  1.7664/(1 +3.0899x-1.8462x2)  , 

where  x  is  defined  as  before. 

The  average  forward  direction  cosine  cos  0  is 
approximated  to  four  significant  figures  on  the 
mterval  0.1  MeV  <  Ey  <  10  MeV  by 

0.660362  -0 .089254/1+0.0985 1 7 p2 

COS  0  ~  - - - 

1.0  -  0.183651/a  +  0.105445/x2 

where 

H  =  tn[o.5+3.7  Ey/(Ey+0.05)]  . 

The  average  recoil  energy  E  is  exactly 
calculable  as 


E  =  Ey ce^!of  , 

where  ae(a)  is  the  Klein-Nishina  Compton  ab¬ 
sorption  cross  section. 


4.  CONCLUSIONS 

In  this  investigation,  an  energy-dependent  elec¬ 
tron  transmission  factor  has  been  used  to  evaluate 
more  accurately  the  Compton  electron  average 
forward  range,  for  incident  gamma  energies  from 
0.3  to  10.0  MeV .  The  results  were  compared  w  ith  a 
popular  previous  approximation  using  an  energy- 
independent  electron  transmission  factor.  The  pre¬ 
vious  approximation  shows  in  comparison  essen¬ 
tially  no  error  for  a  gamma  energy  of  1.5  MeV  , 
7  percent  error  for  0.3  MeV  ,  and  10  percent  error 
for  10.0  MeV. 

For  a  typical  gamma  spectrum  hardened  to  an 
average  energy  of  about  3.5  MeV  by  atmospheric 
attenuation  of  lower  energies,  the  previous  approx¬ 
imation  would  give  Compton  currents  5  percent 
smaller  and  electric  fields  perhaps  10  percent 
smaller  than  for  the  more  accurate  treatment 
presented  here.  Although  the  differences  are  small, 
the  more  accurate  estimate  of  average  forward 
range  presented  here  could  be  incorporated  into 
EMF  environment  prediction  codes  at  virtually  no 
cost  Simple  numerical  fits  to  relevant  parameters 
are  given. 


I 


DISTRIBUTION 


DIRECTOR 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 
1400  WILSON  BLVD 
ARLINGTON,  VA  22209 
ATTN  TIO 

DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION 
ALEXANDRIA ,  VA  22314 
ATTN  DDC-TCA  (12  COPIES) 

COMMANDER 

US  ARMY  RSCH  &  STD  GP  (EUR) 

ATTN  LTC  JAMES  M.  KENNEDY,  JR. 

CHIEF,  PHYSICS  &  MATH  BRANCH 
FPO  NEW  YORK  09510 

COMMANDER 

US  ARMY  MATERIEL  DEVELOPMENT  & 

READINESS  COMMAND 
ATTN  DRXAM-TL  HQ  TECH  LIBRARY 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

COMMANDER 

US  ARMY  ARMAMENT  MATERIEL 
READINESS  COMMAND 

ATTN  DRSAR-LEP-L,  TECHNICAL  LIBRARY 
ROCK  ISLAND,  IL  61299 

COMMANDER 

US  ARMY  MISSILE  &  MUNITIONS 
CENTER  &  SCHOOL 
ATTN  ATSK-CTD-F 
REDSTONE  ARSENAL,  AL  35809 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS 
ACTIVITY 
ATTN  DRXSY-MP 

ABERDEEN  PROVING  GROUND,  MD  21005 
DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN  DRDAR-TSB-S  (STINFO) 

ABERDEEN  PROVING  GOUND,  MD  21005 

TELEDYNE  BROWN  ENGINEERING 
CUMMINGS  RESEARCH  PARK 
ATTN  DR.  MELVIN  L.  PRICE,  MS-44 
HUNTSVILLE,  AL  35807 


ENGINEERING  SOCIETIES  LIBRARY 
345  EAST  47TH  STREET 
ATTN  ACQUISITIONS  DEPARTMENT 
NEW  YORK,  NY  10017 

US  ARMY  ELECTRONICS  TECHNOLOGY 
&  DEVICES  LABORATORY 
ATTN  DELET-DD 
FORT  MONMOUTH,  NJ  07703 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  DC  20305 
ATTN  TITL 

ATTN  RAEV  (G. BAKER) 

ATTN  DDST 

US  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
7500  BACKLICK  ROAD 
BUILDING  2073 
SPRINGFIELD,  VA  22150 
ATTN  DR.  J.  BERBERET 

AIR  FORCE  WEAPONS  LABORATORY,  AFSC 
KIRTLAND  AFB ,  NM  87117 
ATTN  EL,  C.  BAUM 
ATTN  ELT,  W.  PAGE 

DIKEWOOD  INDUSTRIES,  INC. 

1009  BRADBURY  DRIVE,  SE 
ALBUQUERQUE,  NM  87106 
ATTN  C.  JONES 

GENERAL  ELECTRIC  CO. -TEMPO 
816  STATE  STREET 
P  O  DRAWER  QQ 
SANTA  BARBARA,  CA  93102 
ATTN  DAS I AC 

JAYCOR 
P  O  BOX  2008 

SANTA  BARBARA,  CA  93120 
ATTN  W.  RADAS KY 

LAWRENCE  LIVERMORE  LABORATORY 
P  0  BOX  808  ' 

LIVERMORE,  CA  94550 
ATTN  DOC  CON  FOR  TECHNICAL  INFORMATION 
DEPT.  LIBRARY 

ATTN  DOC  CON  FOR  L-94  H  KRUGER 


11 


; 


DISTRIBUTION  (Cont ' d) 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P  O  BOX  1663 
LOS  ALAMOS ,  NM  87545 
ATTN  DOC  CON  FOR  J  MALIK 

MISSION  RESEARCH  CORP . 

P  O  DRAWER  719 
SANTA  BARBARA ,  CA  93102 
ATTN  EMP  GROUP 


R  &  D  ASSOCIATES 
1401  WILSON  BLVD 
SUITE  500 

ARLINGTON,  VA  22209 
ATTN  J.  BOMBARDT 

SCIENCE  APPLICATIONS,  INC. 

P  O  BOX  2351 
LA  JOLLA,  CA  92038 
ATTN  R.  PARKINSON 

US  ARMY  ELECTRONICS  RESEARCH 
&  DEVELOPMENT  COMMAND 
ATTN  TECHNICAL  DIRECTOR,  DRDEL-CT 


HARRY 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


DIAMOND  LABORATORIES 

00100,  COMMANDER/TECH  DIR/TSO 

CHIEF,  DIV  10000 

CHIEF,  DIV  20000 

CHIEF,  DIV  30000 

CHIEF,  DIV  40000 

RECORD  COPY,  81200 

HDL  LIBRARY,  81100  (3  COPIES) 

HDL  LIBRARY,  81100  (WOODBRIDGE) 

TECHNICAL  REPORTS  BRANCH  81300 

CHAIRMAN,  EDITORIAL  COMMITTEE 

CHIEF,  21000 

CHIEF,  21100 

CHIEF,  21200 

CHIEF,  21300 

CHIEF,  21400 

CHIEF,  21500 

CHIEF,  22000 

CHIEF,  22100 

CHIEF,  22300 

CHIEF,  22800 

CHIEF,  22900 

SPOHN,  D. ,  21300 

KERRIS,  K.,  22900 

WIMENITZ,  F.,  20240 

WYATT,  W.T.,  21300  (3  COPIES) 


12 


END 


DATE 

FILMED 


DTIC 


•  L-t'-iJi-;*;.;.  >jiit  iAaiw'MgU fi  • 


